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Simulation of conditional master equations is important to describe systems under continuous 
measurement and for the design of control strategies in quantum systems. For large bosonic sys- 
tems, such as BEC and atom lasers, full quantum field simulations must rely on scalable stochastic 
methods whose convergence time is restricted by the use of representations based on coherent states. 
Here we show that typical measurements on atom-optical systems have a common form that allows 
for an efficient simulation using the number-phase Wigner (NPW) phase-space representation. We 
demonstrate that a stochastic method based on the NPW can converge over an order of magni- 
tude longer and more precisely than its coherent equivalent. This opens the possibility of realistic 
simulations of controlled multi-mode quantum systems. 
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Exciting advances in physics have led to a boom of re- 
search into technologies that exploit fundamental quan- 
tum properties. Such quantum technologies now encom- 
pass more than lasers and superconductors. Indeed, there 
are applications to precision metrology jl] |5] , quantum 
information processing and quantum cryptography [31 H] . 
A key feature of quantum technologies is that they re- 
quire the precise creation, measurement and control of in- 
dividual quantum systems. In particular, measurement- 
based feedback control has shown promise as an effective 
and robust technique for controlling quantum systems. 
The first experiments O Ej and much theoretical work 
[7lfT2] on feedback control of quantum systems have been 
applied to relatively low-dimensional systems. This let- 
ter describes a technique for efficient simulation of large 
bosonic conditional quantum systems that is more than 
an order of magnitude more precise and converges for sig- 
nificantly longer timescales than previous methods, and 
that scales logarithmically with the size of the Hilbert 
space. 

A large bosonic system of particular interest to 
quantum science is the Bose-Einstcin condensate. 
Measurement-based feedback control of BECs and atom 
lasers was first investigated in a single-mode model, 
where a continuous number measurement was used to 
reduce the interaction-induced phase diffusion that lim- 
its single-mode atom laser linewidth [T3]. It was then 
shown that position measurement and feedback on a sin- 
gle trapped atom could bring it to the ground state [9], 
but the proposed measurement scheme was not suitable 
for large atomic clouds such as a condensate. A multi- 
mode quantum field model of a condensate measured by 
an existing experimental technique (phase-contrast imag- 
ing) was then produced, but it could only be solved using 
a semiclassical approximation [Til IE] ■ Analysis of the 
linewidth of a multimode atom laser undergoing feedback 
will require a viable stochastic method for conditional 



quantum states that can deal with both high nonlinear- 
ities and number-like measurements. The number-phase 
Wigner function method fulfills both of these require- 
ments. 

The most effective methods for dynamic simulation of 
high-dimensional bosonic quantum systems are stochas- 
tic techniques based on phase-space representations |14l 
[15] . Each stochastic method is derived from a specific 
phase-space representation, which is akin to the choice 
of a basis for the Hilbert space. Naively, these tech- 
niques require memory and computational resources that 
scale logarithmically with the size of the Hilbert space. 
Practically, the overall computational efficiency is sys- 
tem dependent, and strongly depends on how well the 
underlying phase-space representation matches the nat- 
ural basis for the quantum system under consideration. 
The most commonly used stochastic simulation meth- 
ods are based on phase-space representations that use 
Gaussian states. These methods have enabled the sim- 
ulation of quantum optical [TSl US] , atomic [131 [TT] [TB] , 
and fermionic quantum fields |19| . In particular, stochas- 
tic methods have been used extensively in the field of 
quantum-atom optics, where dilute atomic gases can 
be cooled to produce BECs and atom lasers [20H22] . 
The two most successful varieties are based on posi- 
tive P (P"*") and truncated Wigner (TW) representa- 
tions. P"*" is an exact technique, but requires a dou- 
bling of the phase space that often leads to instabili- 
ties |23j . Truncated Wigner is an approximate technique 
that typically has significantly longer convergence times 
than P"*". However it makes an uncontrolled approxima- 
tion [21], and may therefore converge to incorrect so- 
lutions. Both of these methods, along with all other 
coherent-state based representations, experience difficul- 
ties dealing with large number-conserving nonlinearities, 
as the underlying Gaussian basis becomes inappropriate. 
Such large number-conserving nonlinearities are typically 
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the dominant energies in confined cold atomic systems 
by a couple of orders of magnitude. Recently, we in- 
troduced a new stochastic method based on a number- 
phase Wigner (NPW) representation [5S], that provides 
a non-approximate method for simulating large number- 
conserving nonlinearities. It was found that this dramat- 
ically improved the convergence of simulations of these 
highly nonlinear systems. 

Modelling highly nonlinear systems undergoing contin- 
uous monitoring and feedback requires the simulation of 
a conditional quantum state. In a recent paper [2B], we 
extended stochastic simulation techniques to apply to a 
class of conditional quantum systems. Continuous mea- 
surement of a quantum system can have a dramatic effect 
on its dynamics. In fact the choice of measurement can 
even be used as a controlling mechanism by itself [27] . It 
is therefore unsurprising that the appropriate choice of 
phase-space representation is heavily influenced by the 
choice of measurement, as it may drive the conditioned 
system towards a state that is simpler to describe in a 
particular representation. Also, not all measurements 
automatically produce a usable method. Our previous 
paper demonstrated how to unravel a particular form of 
stochastic Fokker-Planck equation (SFPE) [26] , and it is 
only possible to generate SFPEs of this form with par- 
ticular combinations of measurement schemes and phase- 
space representations. In particular, methods based on 
coherent state representations are badly suited to mea- 
surements involving number-like observables rather than 
quadrature-like observables, which are prevalent in atom 
optics. In this letter, we show that the NPW representa- 
tion produces dramatically superior results than coherent 
state based representations for these calculations. 

A common quantum atom-optical system under mon- 
itoring is governed by the conditional master equation 

dp^~i[H,p] dt + ^V[L,]p dt + J2'H[^P dW^. (1) 

i i 

where dW is an Ito Wiener increment; I?[c]p — cpc^ — 
\{c^cp + pc^c); y.[c] = cp + p& — Tr[cp + pc^]p] H is 
the Hamiltonian and contains the contributions from ki- 
netic, potential and many-body interaction energies; and 
Li = J dxtp^ {x)Li{x)i[:{pc) is the measurement operator 
where Li(x) € are the measured density moments 
of the multimode object. The only restriction we have 
applied to the measurement operators Li is the order 
of the field operators which we note are 'number like' 
algebraically. This form is the lowest order number- 
conserving interaction possible for a multimode system. 
Measurements that are lower order with respect to the 
field operators will not conserve number and may be 
suited to traditional coherent-state representations, but 
in many cases these systems may be treated using an- 
alytic techniques like the Kalman filter, making simula- 
tion less important. Number-conserving measurements 
are quite common in engineered monitoring of BECs 



[TTl [T^ [25H55] , as any phase-sensitive measurement re- 
quires the existence of an atomic local oscillator to use 
as a phase standard. Thus, the efficient simulation of 
Eq. ([T]) will be relevant to a wide variety of atom-optic 
systems, including all those involving current experimen- 
tal detection schemes. 

To compare the performance of the NPW represen- 
tation to coherent methods when simulating conditional 
master equations of the form Eq. ([T]) we require a veri- 
fiable solution for comparison. Unfortunately, Gaussian 
analytic techniques commonly applied in quantum con- 
trol are not appropriate for Eq. ([T]). This can be under- 
stood by noting the measurement operator Li is second 
order with respect to the field operators, which generates 
non-quadratic terms. Thus Gaussian analytic techniques 
such as the Kalman filter are not guaranteed to be exact 
[34] . and we are forced to integrate the master equation 
directly to generate a benchmark for comparison. This 
restricts us to looking at single mode systems, as direct 
integration is not scalable to multimode systems. The 
single mode problem that is algebraically equivalent to 
the multimode Eq. ([I]) is 

dp = 7X> [a"^ d]pdt + -fC [a^ d]p dt + y/^H [a'f d] o dW. (2) 

where odW is a Stratonovich 

Wiener increment and C[c]p ~ 
-~\ {(?p + 2cpc^ + p{c^f - Ti[(?p + 2c0 + p(ct)2]p) + 
(cp + pc' — Tr[c/5 + pc^]pj Ti-[cp+pc^] is the Stratonovich 
correction superoperator. This master equation is of a 
system undergoing continuous coUopse under a number 
measurement. 

The scalability of stochastic techniques for solving con- 
ditional quantum dynamics has already been demon- 
strated in ,26l , but we aim to investigate the effect 
of choosing different representations. We use master 
equation Eq. ^ to compare the performance of leading 
coherent-based scalable stochastic methods to the NPW 
representation. The convergence of these techniques is 
compared to a direct integration of the master equation. 

We start our analysis with the coherent state based 
representations P+ and TW. The success of these tech- 
niques have been primarily concerned with Hamiltonian 
and decoherence evolution of BEG and quantum-optical 
systems. Starting with P"^, we now investigate the ap- 
plicability of these techniques on a conditional master 
equation. Using the correspondences in |15j we can con- 
vert the master equation ^ to 

dV{a) = {7 [daa (1 + 2|a|2 - 2Ep [\a\^]) 

+dc,,a* (l + 2|ap-2Ep 

_2(|a|2 + |a|4-Ep [jap] - [|a|4] ) 
+4Ep [lap] (|ap-Ep [jap])] dt 
+y/l[-daa - da' a* 

+2(|a|2-Ep [lap])] odI^}7'(a), (3) 
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where Vi^a) is the P-representation quasi-probability 
distribution which reproduces normally ordered mo- 
ments of the master equation Here Eq [/(x)] = 
/ c?x'/(x')Q(x') is our notation for taking the expecta- 
tion values of a function /(x) with respect to the quasi- 
probability distribution Q(x).We immediately note that 
this equation contains non-positive definite diffusion that 
must be simulated by doubling the phase space. Thus P+ 
techniques are required. This representation can then be 
unravelled into the following set of stochastic equations: 



(4) 



da = -27a {/3a -Ef [(3a]) dt 

+ ° + *^^2 + dW); 

d(3 = -27/3 {Pa-Ef [f3a])dt 

+ Vll^ ° (--i-dVi + idV2 + dW); 

duj = -27W {^a + (3^a^ - 2(3aEf [(3a]) 
+ 2y/^Lol3a o dW. 

Where dVi and dV2 are the set of 'fictitious noises' 
that are averaged over to obtain the weighted averages 
Ef [/(x)] = ^ii^if{^i)/J2i^i- ^'^^ more details on the 
techniques used to unravel equation Q into Q and how 
to simulate them see [55]. Eqs. Q will be used to bench- 
mark the unravelling of Eq. (pi using coherent-based 
methods. 

We continue our analysis with the TW representa- 
tion. Using the operator correspondences given in [15], 
we can write the master equation for the Wigner quasi- 
probability distribution yV{a.) as 



dW{a) = \ 7 



-Idiiar 



a 



Eyy 

^ |4 



2(|a|4-Ew [|«|^]) 
4Ew [lap] {\a\^~Ew [|ap])] dt 



+ 



1 



2(|ap-Ew [lap]) 



odW^yV{a) 



Finally we consider the number-phase Wigner repre- 
sentation. The NPW was first derived in [25] and was 
used in the simulation of large nonlinear equations. We 
now consider its applicability for use on conditioned large 
atom-optic systems. Using the operator correspondences 
given in |25| we get the following equation 



(5) 



Note that the first term contains higher order derivatives 
and a truncation is required in order to obtain a stochas- 
tic unravelling of Eq. ([5| . Note also that a P"*" style ex- 
tension of the phase space would be required to simulate 
the diffusion in the conditioning term. Traditionally the 
Wigner representation is guaranteed to produce strictly 
positive-definite diffusion [15], but this is under the as- 
sumption that the calculus increment is positive as is the 
case with dt. Unfortunately this assumption does not 
hold with the dW increment. A new 'positive' Wigner 
representation could be derived by analogy to the P+ rep- 
resentation, but the higher order terms would still need to 
be truncated. This would make this hypothetical repre- 
sentation both approximate and doubled in phase space, 
which would make it unlikely to compete with P"*". Thus 
it is not worthy of further investigation. 



dM{n, = 7 



dl-2{n^-EM[n^]) 



+AEj^[n]{n-Ej^[n]) dt 

+2^[n-E^r[n]]odW^^^{n,<j>), (6) 

where Af{n, (j)) is the NPW representation that produces 
a complete set of moments of the master equation as 
outlined in [25j. We next unravel Eq. Q using [26] to 

dn — 0; 

dfj) = ^/ldVl] 

duj = 7w(-2n^ -I- 4E/ [n] n)dt + ^bjn o dW. (7) 

Note we did not need to apply any truncations or dou- 
ble the phase space. The simplicity of the equations ^ 
compared to Q show how an appropriate choice of rep- 
resentation, NPW in this case, can greatly reduce the 
complexity of the evolution, just as an appropriate choice 
of basis can simplify analysis of other quantum problems. 
We can now compare the performance of the NPW rep- 
resentation to P+ by integrating equations Q and ([7|, 
respectively. 

The simulations comparing the P+, NPW and a direct 
integration of the master equation are shown in Fig. [Tj 
The number-phase Wigner representation converges for 
the longest time interval. In fact it converges until a com- 
plete collapse into the correct number state. As this is 
the steady state of the equation we expect the number- 
phase representation to converge indefinitely. Not only 
is the NPW more accurate it is also significantly more 
precise. The increase in precision in turn improves the 
accuracy of the evolution, as simulation of the condi- 
tional master equation uses an estimate of the observ- 
able Tr[c/) + pc^], thus lack of precision results in a lack 
of accuracy in the long term. This dynamic instability 
is not seen in non-conditional master equation evolution, 
and makes the precision of stochastic techniques consid- 
erably more important in these problems. 

The results show that the NPW-based simulations 
are significantly better than the simulations based on 
coherent-state representations for the conditional master 
equation described in Eq. ([2]). The simulations promise 
to be stable enough to consider long-term behaviour of 
systems and model the effects of feedback strategies. This 
high level of convergence can be explained by noting the 
suitability of the basis underlying the representation to 
the measurement eigenstates. Importantly, due to prac- 
tical difficulties in producing stable atomic local oscil- 
lators, any current detection scheme used in quantum 
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FIG. 1. (Color online) (a) Plot of number versus time (in 
dimensionless units) for a single measurement run, compar- 
ing the direct integration of the master equation to the NPW 
and methods. Simulations had an initial condition of a co- 
herent state with an amplitude of 10. Numerical integration 
of the master equation (|2| is plotted with a solid green line, 
the NPW ^ is plotted with a dashed red line and P+ Q 
is plotted with a dot-dashed blue line. Uncertainty in both 
stochastic methods is plotted with dotted lines. The P^ be- 
comes divergent around t = 0.15, and is not plotted beyond 
this point. A close-up view of this divergence is shown in the 
inset. Part (b) shows the accuracy, defined as the difference 
between each stochastic result and the master equation so- 
lution, and the precision, defined as the standard deviation 
of the averages, for each stochastic method on a logarithmic 
scale. The NPW accuracy is described with a solid red line, 
the NPW precision is a dashed red line, the P^ accuracy is a 
dash-dotted blue line, and the P"*" precision is a dotted blue 
line. The NPW is considerably more precise and is convergent 
for at least an order of magnitude longer than the competing 
P^ representation. The numerical integration was performed 
by using the open source software package XMDS |35) . 



gases is of a form that is suited to the NPW method. 
The NPW method described in this letter is also the 
only simulation tool that is deterministic for the strong 
number-conserving nonlinearities that are present in such 
systems, and is therefore the only suitable candidate for 
simulating conditional states of ultracold atomic gases 
for feedback or state estimation. 
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